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We explore the consequences of the existence of a very large number of light scalar degrees of 
freedom in the early universe. We distinguish between participator and spectator fields. The former 
have a small mass, and can contribute to the inflationary dynamics; the latter are either strictly 
massless or have a negligible VEV. In N-flation and generic assisted inflation scenarios, inflation 
is a co-operative phenomenon driven by N participator fields, none of which could drive inflation 
on its own. We review upper bounds on N, as a function of the inflationary Hubble scale H. We 
then consider stochastic and eternal inflation in models with N participator fields showing that 
individual fields may evolve stochastically while the whole ensemble behaves deterministically, and 
that a wide range of eternal inflationary scenarios are possible in this regime. We then compute one- 
loop quantum corrections to the inflationary power spectrum. These are largest with N spectator 
fields and a single participator field, and the resulting bound on N is always weaker than those 
obtained in other ways. We find that loop corrections to the N-flation power spectrum do not scale 
with TV, and thus place no upper bound on the number of participator fields. This result also implies 
that, at least to leading order, the theory behaves like a composite single scalar field. In order to 
perform this calculation, we address a number of issues associated with loop calculations in the 
Schwinger-Keldysh "in-in" formalism. 



Many candidate theories of fundamental physics predict the existence of large numbers of scalar degrees of freedom. 
Classically, if these modes are not excited they play no role in the cosmological dynamics. Quantum mechanically, 
however, light scalar modes fluctuate with an amplitude set by the Hubble scale H and, since everything couples 
to the graviton, they contribute to loop corrections. Thus, while adding N light scalar modes need not change the 
classical dynamics of the early universe, we expect quantum contributions that scale with N. In addition, if these 
fields are in thermal equilibrium with the rest of the universe their contribution to the effective number of degrees of 
freedom modifies the relationship between density and temperature. For a given H we can therefore find an upper 
limit on N, above which these modes dominate the cosmological evolution. 

Consider a scenario with N scalar fields, 



where the total potential V = ^/Vj (</>/) is a sum of N uncoupled potential terms. 1 A field 0/ is "light" if 
drVi/dcpr/2 = rrij <C H 2 . We make a distinction between participator and spectator fields. The latter are either 
massless because Vt is strictly zero, or their vacuum expectation value (VEV) is very small. Conversely, participator 
fields have small but non-zero mass and sufficient VEV for them to help drive inflation via their contribution to the 
overall density. 

Perhaps the simplest route to a meaningful bound on N is to note that all these fields undergo quantum mechanical 
fluctuations of order 6(f>i ~ H/2ir. During inflation these fluctuations freeze out as classical perturbations at scales 
larger than the Hubble length, 1/H. Each field has gradient energy (V0) 2 /2, which counts towards the overall energy 
density. The gradient energy thus scales like N(5(f>/Sx) 2 jl ~ NH 4 /8ir 2 . Given H, the energy density is provided by 
the — Einstein equation, H 2 — p/3M 2 . For self-consistency, the gradient contribution must be much smaller than 
other contributions to p, so 
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1 We use the reduced Planck mass M£ = 1/SnG throughout. 
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If inflation occurs at the GUT scale, then M p /H ~ 10 5 and N <§C 10 10 . This bound can be derived by many routes 
(e.g. [3, 0), and appears to be robust. 2 

One can also consider bounds on N from loop corrections to the gravitational constant. Since gravity (and thus 
the graviton) couples to all fields, matter loop corrections can renormalize its value @, Q- Veneziano [7|] argued that 
in order for the effective value of Newton's constant to be greater than zero we need 

Ml 

(3) 

where A is the scale of the invariant UV cut-off, for example the mass scale of the N stable fields. Likewise, Veneziano 
Q points out that this bound prevents potential violations of the holographic bound on the entropy density that can 
be encountered when the total number of species grows without limit [8, 9], the so-called "Species Problem". More 
recently, Dvali noted that in the presence of a large number of species, the Veneziano mechanism weakens the 
gravitational coupling by a factor of l/N. Setting N ~ 10 32 and working at the TeV scale to satisfy the bound of 
equation (|3|) solves the hierarchy problem, provided some ultraviolet completion of the standard mode can produce 
the requisite value of N. Dvali also provides an alternate non-perturbative derivation of equation based on the 
consistency of black hole physics. 

In Section |TT] we begin by considering stochastic inflation jTl], [T^, EH with multiple degrees of freedom. Many 
simple inflationary potentials possess a range of field values for which the potential is safely sub-Planckian, while 
the stochastic motion of the field dominates the semi-classical rolling. With a single field, a stochastic phase is 
necessarily eternal, since the inflationary domains perpetually reproduce themselves. However, once N > I, more 
complicated scenarios become possible. Specifically, with N participator fields we can implement assisted eternal 
inflation without any single field acquiring a super-Planckian VEV. This scenario is a variety of "assisted inflation" 
where the inflaton is a composite of many individual fields (l4| . Secondly, we find solutions where individual fields 
move stochastically, but a well-defined composite field rolls smoothly towards its minimum. Finally, we find models 
where a field (or fields) evolves semi-classically, along with other fields that move stochastically. If these fields have 
symmetry breaking potentials they yield an apparent "multiverse" of disconnected bubbles. However, the stochastic 
phase ends globally, as the semi-classically evolving field gives a natural cut-off to what would otherwise be an eternally 
inflatin g un iverse , providing a potential toy model for studying the well-known measure problem in eternal inflation 

[3 EE EE E3, El, [H, . 

All fields couple to the inflaton gravitationally, and thus contribute loop corrections to the inflation potential, via a 
coupling of order (H/M p ) 2 . This is necessarily a small number during the last 60 e-folds of inflation, given that H fixes 
the scale of tensor fluctuations and is bounded by the absence of an observed B-mode in the microwave background. 
Consequently, any single loop makes a tiny correction to the inflaton propagator. However, this contribution is 
amplified by V, the number of species that can flow round the loops, and in Section IIIII we compute the relevant 
one-loop corrections to the inflaton propagator. We consider two limits - N spectator fields with a single inflaton, 
and the N- nation case, with N participator fields [U • With N spectator fields we find an upper bound on N 
similar to those of 0, [23| • We describe these results in Section IIIII relegating many details to Appendix [SJ We work 
in the "in-in" formalism [24], [25|, HH, H3 j which has been turned into an extreme ly p owerful tool for studying higher 
order corrections to cosmological correlations by Maldacena [28| and Weinberg [29| ]. This calculation requires us to 
develop the fourth order interaction Hamiltonian for a theory of inflation with N uncoupled scalar fields, which we 
present in Appendix [B] and which will have applications beyond the current calculation. Moreover, it turns out that 
there are some subtleties with the computation of loop corrections in the in-in formalism which we also clarify in the 
Appendices. Finally we conclude in Section [TVl 



II. STOCHASTIC V-FLATION 



In simple versions of V-flation, one has N identical participator fields, and inflation emerges as a cooperative 
phenomenon. We assume that the overall potential is the sum of the individual fields' potential terms, and that 
cross terms are absent. Each field only feels the "slope" of its own potential, but the corresponding friction term in 
the field's equation of motion is still proportional to H, and thus grows with N. Interestingly, a similar scaling also 
applies to stochastic inflation, which occurs when the inflaton evolution is dominated by quantum fluctuations, rather 
than semi-classical rolling |lf|, |l2j, ll3| . The individual fields have fluctuations of order H, so this amplitude will grow 
relative to the semi-classical evolution as N is increased. 



2 in g,a this limit on N is derived in the context of N-flation but in reality it applies to any scenario with N light fields. 
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Let us begin by looking at N fields with identical potential terms and initial VEVs. Assuming slow roll, 

where = Vr (4>i) is the potential for any one of the fields. The amplitude of the quantum fluctuations of the I— th 
field is 

TJ 

S(/>I,q = <50 g = — , (5) 

where the first equality reflects our assumption of identical fields. From equation ((4]) 



(6) 
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Conversely, the distance travelled by 4>i m a single Hubble time is 

.1 VI 1 M;?V7 
H 3H H NV 

where V[ = dV/d(f>i — dVi/d(f>i. The defining condition for the stochastic inflation is that 8<j> c < 6<f> q Forming 
the ratio of these terms gives 

/ - \ 3/2 



12tt 2 \M 2 ) V/' 

and stochastic inflation therefore occurs whenever the above expression is larger than unity. This ratio increases with 
N if everything else is held fixed. This is to be expected, since boosting N lifts the overall energy density, and thus 
the amplitude of the quantum fluctuations. Conversely, the semi-classical rolling slows with N, as we are increasing 
H while holding V[ constant. 

The above discussion holds for a generic potential, but now consider N quadratic potentials with identical masses, 
Vi = m 2 cf)j/2. Assisted inflation can be described in terms of an effective single field, ip which, for quadratic potentials 
is 0/ = ip/VN and [HI HI 



N N<j)'j m _ I N <p 2 to ^ 
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Interestingly, this expression retains an explicit dependence on TV, whereas N-flation ends at a fixed value of ip, 
independently of N. Before proceeding, let us consider some specific numbers. The principal virtue of N-flation is 
that we can build a GUT-scale inflation scenario without invoking trans-Planckian VEVs. Consequently, if <pj ~ M p 
the critical value of V at which the fields can move stochastically is 

V = (967T 2 ) 1 / 3 ( ^ « 10 ( V3 . (10) 
\ m J \ to / 

The ratio m/M p is fixed via the amplitude of the perturbation spectrum, and is around 2x 10~ 4 |2lll2^ . Consequently, 
we might see the onset of stochastic motion in the fields if N ~ O(10 4 ). This is an order of magnitude larger than the 
number of fields one expects in N-flation, on the basis of the likely number of two-cycles (from which the N axions 
are derived) one can find in a realistic Calabi-Yau [HI, [HJ , but is much less than the absolute upper limit on V given 
by equation ([2|). To satisfy this requirement, we need H 2 fa Nm 2 /6 (assuming <f>i — M p ), in which case 

iV<V6^~10 5 . (11) 

TO 

The theoretical upper limit on N is roughly an order of magnitude larger than the value required for the individual 
fields to move stochastically. 3 Conversely, slow roll ends when tp ss \/2M p or tfi ~ y/2/NM p . For self-consistency, we 



3 Huang et. al 
introduced in [31 



30] has argued that choatic eternal inflation with large N fields is ruled out by the so-called "weak gravity conjecture" 
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expect the N fields to be moving semi-classically at this point, so we can derive a very weak upper bound on N by 
setting ip « y/2M p in ©, namely N < 24ir 2 M 2 /m 2 ~ 10 10 , far above even the weak bound of equation Since the 

last sixty e- folds occur at values of (p a few times larger than y/2M p , we need not worry that the fields are moving 
stochastically over astrophysically interesting scales. 

With a single field, the onset of stochastic inflation is synonymous with the amplitude of the density fluctuations 
exceeding unity [32] . For iV-flation the density fluctuations have an amplitude [2l|, [22j 



o 1 / 2 1 tp 2 m 

K " = \ 967T 2 ~M 2 ~Mp~ ' (12) 

which is y/N smaller than the critical ratio for the onset of stochastic motion by the individual fields, so at the 
threshold for stochastic motion, P^ 2 ~ 1/ \/N <C 1. In this case stochastic inflation is not synonymous with eternal 
inflation. As described above, in one Hubble time each field makes a random jump of ~ 8<p q with undergoing 
semiclasical evolution 6<f> c . For a given field, the sign of 5<p q is random, while the semiclassical motion always points 
in the downhill direction. The average field thus has a stochastic fluctuation 8<j) q /yJ~N , since this is the mean of N 
signed, random variables. Consequently, the collective motion of the ensemble of fields remains deterministic, and the 
density of the universe will decrease monotonically with time unless 8<p q > y/N5(f) c . In this case the density fluctuation 
is boosted to order unity, and the average density can increase inside a given Hubble patch, so inflation is not just 
stochastic but eternal. Since H oc p 1 ^ 2 , we are in the intermediate range where inflation is stochastic but not eternal, 
5<j) q /5<j) c diminishes with time. In this case the stochastic motion of the individual fields will eventually become 
subdominant and for any reasonable value of N the inflationary dynamics will be well-described by the semi-classical 
motion alone. 

From a practical perspective, a period of stochastic inflation in a simple multi-field model need not modify the 
observable properties of the universe. Certainly in the case described above, the stochastic motion would cease well 
before P-^ 2 ~ 10 -5 unless N is very large. However, since the stochastic motion necessarily increases the variance 
in the individual field values this phase may have an impact on the likely initial spread in the values the which 
can have an impact on the inflationary observables in iV-flation. However, there is no clear expectation for the likely 
initial values of the (pi and without this the impact of any stochastic evolution cannot be evaluated. 

On the other hand, if the individual potential terms do not all have a single well-defined minimum, any phase 
where one of more of the individual fields moves stochastically could have a substantial impact on the inflationary 
phenomenology. For instance, in the case of iV-flation, the N fields are actually axions and thus have periodic 
potentials. The individual m 2 (f> 2 terms arise from assuming that each field is close to it minimum, when measured 
relative to the scale on which the underlying potential is periodic, and then Taylor expanding. However, if we 
imagine an initial state when the N fields (or even some subset of them) are close to the maxima of their cosine 
potentials two effects will occur. The first is that these fields will have a small 8<p c , since the corresponding V[ will 
be very small in the vicinity of an extremum, making it more likely that these fields move stochastically even if 
other fields are dominated by the semiclassical rolling. Secondly, with fields moving stochastically near the maxima 
of these potentials, in individual Hubble domains in which these fields do evolve away from the peaks the fields will 
then roll towards different minima. In eternal inflation each patch in which the stochastic motion ceases would be 
identified as a "pocket universe" and this scenario thus produces pockets with many different vacua, depending on the 
symmetry breaking pattern of the individual axions. However, if enough fields are rolling semi-classically, H is strictly 
decreasing, and inflation is not future-eternal. In this case we can end inflation globally, and count the number and 
type of distinct domains created during the stochastic phase. Consequently we now have a toy model that initially 
resembles an eternally inflating universe, but in which there is a natural late-time cutoff which removes the infinities 
which otherwise can prevent one calculating the relative creation rates for different types of pocket. This system may 
thus prove useful for testing different "measures" for eternal inflation [IE EE EE El El 

EE 113]' an( i we wm examine 

this issue in a separate publication. Conversely, once inflation has ended, these different domains will eventually 
merge and the late-time universe (in the absence of a cosmological constant term) will be composed of domains with 
different vacua, separated by domain walls. 



III. ONE-LOOP QUANTUM CORRECTIONS 



In the previous section, we considered stochastic field evolution in cosmological scenarios with many degrees of 
freedom. We now assume the field evolution is well-described by the semi-classical equations of motion, and compute 
loop corrections to the perturbation spectrum. Scalar loop corrections to the inflationary spectrum in single field 
models have been widely discussed [29ll33l. [34l l35l. l36l. \3l\ . Likewise, the correction from graviton loops was computed 



5 



in [38[ . 4 Weinberg [2^| looks at loop corrections to a single inflaton in the presence of many massless spectator fields, 
and we now generalize this result to models with N participator fields. For clarity, we present our results in this 
section, and discuss technical aspects of the calculation in Appendix |A] 

We previously made the distinction between single field inflation with N spectator fields, and assisted scenarios with 
N participator fields, of which A-flation is the most interesting example. The field dynamics differs between these 
cases, since only the participator fields have non-zero VEV and contribute to the vacuum energy that drives inflation. 
However, both classes of field can contribute loop corrections to the 2-point function or perturbation spectrum - and 
we will see that the forms of these contributions are different. The immediate concern is that loop corrections result in 
a new "species problem" , as pointed out in 0, HJ, [39| . Generically, if we insist that quantum corrections to the Planck 
scale are small, we need ~ A(A[/y /M p ) 2 to be small, where A^y is some UV cut-off scale. Having computed the loop 
corrections to the 2-point correlation function, we can then determine whether these also provide an constraint on N, 
as a function of Aj/y /M p . Surprisingly, with N participator fields, there is no bound on A. On the other hand, with 
N spectator fields, we obtain slightly weaker bounds on N as compared to the standard arguments. In what follows 
we look at the two limiting cases, firstly A-spectator fields with a single inflaton and then N participator fields. One 
could easily generalize our result to the case where one had a mix of both spectator and participator fields. 



A. One-loop corrections of N Spectator fields 

Consider inflation driven by a single field field 0, with N spectator fields, o~i where / runs from 1 to A, 

(13) 



f M 2 1 1 

S = J dx^g -fR~^ (dcpf + V{4>) - 2 (^) 



In the discussion below, repeated indices are not summed over unless explicitly specified. Unlike the A-flation case, 
the spectator fields remain invariant under the shift cj — > 07 + <5j. Any initial kinetic energy possessed by these fields 
decays away rapidly, since pi cx (<7j) oc a~ 6 . Thus the only contribution to the background energy density is the 
inflaton potential, V{cf>). This situation has been discussed by Weinberg 29], who found that the one-loop two- vertex 
quantum corrections modify the power spectrum equation (|24j) by a term of order (H/M p ) 4 \nk per field. Thus, for 
N ai fields the first order correction to the power spectrum is 

fO- 1 N lHl e l nk (141 
k 4(2tt) 3 6M* 6 *' [ ' 

The one-loop corrected power spectrum is then 

1 1 H 2 ( , 7T H 2 \ 

^w^( 1 + ( ci+ V^ ? lnfc )< (15) 

where C\ = —2tt/3 is the one-vertex self-correction term which we compute in the appendices. What about the one- 
loop one- vertex loops of the a fields around the inflaton? As we explained in detail in Appendix [X] and later in Section 
IIIIB1 there is no non-scale-free <jj one-vertex correction to the <fi propagator, and hence Weinberg's computation is 
complete modulo the inflaton self-correction. 

Requiring that the one-loop corrections do not dominate the "tree-level" two-point correlation, we obtain the bound 

Ml 1 

N<j£-. (16) 

For successful slow roll inflation e <C 1, so this bound is necessarily weaker than that obtained from gradient energy 
considerations. Interestingly, the "tree-level" power spectrum for the tensor modes for single scalar field inflation is 

2 



4 In this work we are primarily interested with the scaling behavior of the loops as one increases the number of fields. Since there are 
only 2 graviton modes, loops involving gravitons cannot scale with the number of fields N (even though individually they are of the 
same order), and hence we neglect them. 

5 The numerical factor in [29j differs from that of equation (|I4[I , As we explain in Appendix |A] there is an extra contribution, related to 
the contour in the time integral needed to pick up the "in" vacuum. We thank Steven Weinberg for a useful discussion of this point. 



FIG. 1: Feynman diagrams for the three- and four-point one-loop correction for the / field correlator. With TV participator 
fields, the first two terms apparently scale as TV, due to the summation of J, but each loop is suppressed by e/TV, arising from 
the extra coupling. At leading order in slow roll, the third term only has contributions from the self-interaction I = J term, so 
the corrections do not blow up as TV becomes large. Here we have denoted I fields by a solid lines while dotted lines represent 
J fields. 

Hence any measurement of P gw would effectively put an observational upper bound on TV, 

TV < 1/P gw . (18) 



B. One-Loop corrections with TV Participator Fields 



We now turn to the case with TV participator fields, N-flation, which is based on the action of equation ([T]), and we 
perturb each individual field, 



Qi 



(19) 



where <f>i is the homogeneous background solution and Qi is the perturbation. In the discussion below, upper case 
Roman letters I, J, ... label the fields, and we have a flat target space X 1 = Xj with the summation convention 



X I Y I =J2X I Y I . 



(20) 



The power spectrum produced by TV uncoupled inflating fields is [4C 



1 
TV 2 " 



E(r) (QiQi) 



(21) 



The 1/TV 2 factor in front reminds us that the total power spectrum is not simply a sum of the individual power spectra 
of the fields; the power spectrum is itself the expectation value of a set of random variables. In general, 



(QiQi) = (QiQi)o + (Q 



iQih + (QiQih..., (22) 

where the numerical subscript denotes the order of the expansion, or the number of vertices. The uncorrected power 
spectrum is simply 

1 



leading to the primordial power spectrum, 



(QiQi)o 



Pi 



2(2vr) 3 



H 2 



4(2tt) 3 M 2 TVe/ ' 



(23) 



(24) 



Our goal is to compute the one-loop corrections to the power spectrum for each field (QiQi}i—i 00 p generated by all 
TV fields. To compute the one-loop correction, we use (|A1[) . with the appropriate interaction Hamiltonian H- lat for 
each field. 6 The three-point interaction Hamiltonian is, to leading order in slow roll [33l |. 



dfx 



3 3 



i.j 



i.j 



(25) 



All the fields below are the taken to be Hcisenberg fields unless otherwise noted. 
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Application of equation (|A4[) leads to three sets of diagrams generated by the two interaction terms and their 
cross- interaction, the 1st two diagrams of Fig. [1] where each vertex corresponds to an interaction via one of the terms 
above. We also need the four-point action, equation (|B6|) derived in Appendix IBl 



d xa 



^ Y d l Qjd l Qjd- 2 (d J Q I d,Q I + Qid 2 Q!) 



4Ha? 



i.j 



\r Y QjQjd~ 2 &QAQi + Qid 2 Qi) 



4H 



j. i 



+-Y d ~ 2 ( d jQj 9 jQj + Qjd 2 Qj)d~ 2 {djQidjQi + Qid 2 Q!) 



i, J 



(26) 



where 0% is given in equation ()B7|) and repeated lower roman indices ■ ■ ■ are summed using the Euclidean metric; 
a i°i = X)i j 3ij a ibj- At first order this interaction generates the third diagram of Fig. [T]via Eq. (|A3|) . For simplicity, 
we ignore the self-interaction term of the two- vertex term, as we expect that its correction to be of the same order as 
the other N — 1 corrections jU, HU . 

Note that in deriving the above expression we have made use of the substitution Hint = — Ant) where Ant is the 
Lagrangian density computed to the appropriate order in the perturbation. This substitution is trivial with non- 
derivative interactions, but it is more subtle in the presence of derivative interactions. However, in Appendix [Cl we 
show that we are effectively ignoring a correction that is at most of order 0(e), 



0(e), 



which can be discarded to leading order in slow-roll. 

Each of the field Q 1 can be expanded in Fourier modes and their respective creation/anihilation operators, 

Q/(x,i)= / ^ke A -'\ ai 0tM(t) + 4(--k)U^(t) 



(27) 



(28) 



where the functions U£(t) are solutions of the equations of motion obtain from varying the second order action, Eq. 
(|B4I) . with respect to the field Q 1 . The a 7 (k) satisfy the usual commutation relations, 



[a 7 (k),4(k')]=*(k-k , )#. 



(29) 



With these ingredients, we make the final simplifying assumption that the fields are identical, which means that we 
simply need to compute the one-loop correction from one of the J terms in equation (|25p and then multiply them by 
N — 1 to get the total correction per field /. We leave the details of this computation to Appendix [A"l 

It turns out that the only the two-vertex diagrams and the one-vertex self-interaction diagram contribute the 
physical logs while non- self-interaction diagrams of the one- vertex loops contribute polynomial ultra-violet divergences 
which we assume are absorbed by renormalization. This is a consequence of the symmetry of the original action where 
the potentials are uncoupled and the kinetic terms are SO(N) symmetric. This can be understood as follows. The 
diagrams in Fig. (J2I) have metric "graviton" propagators and loops hidden inside them - the higher order interactions 
that appear perturbatively are mediated by gravity. Howe ver, we have chosen a gauge where the metric perturbation 
vanishes. If we choose a different gauge (Newtonian gauge [43| for example) these interactions are manifest, and each 
one-vertex diagram receives contribution from terms like 



(30) 
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where the wiggly lines represent the graviton lines, the solid lines are the field perturbation lines and the ellipses 
denote terms that are higher order in metric perturbation or the slow roll parameters. The first "balloon" diagram 
diverges polynomially; the loop that runs around itself does not possess a scale since it is not directly connected to 
an external leg. Both vertices of the second term must obey the symmetry of the original action, which forces I = J. 
It is clear from this perspective why only the self interaction one- vertex loops can contribute to physical obervables. 7 
After regularizing the results and assuming that we can always find suitable counterterms to cancel the gravitational 
backreaction and UV divergences from loops independent of external momenta, the correction of J loops to the / 
propagator is given by equations (|A16[) and (|A29|) . The total correction is thus 



(<9/<9/)i-loop — 



1 H 4 

2(2tt) 3 



ci 



M 2 H 2 



In k. 



(31) 



where c\ = —2ir/3 and c 2 — (2017/240)7r, which arise from the one-vertex self-interaction and the two-vertex loops 
respectively. We have set N — 1 — > N, since N ^> 1, so the two-vertex self-interaction is not important. During 
7V-flation, each participator field is corrected by its N — 1 counterparts in the two-vertex loop, and by itself in the 
one- vertex loop. However, in equation (13ip . the two- vertex corrections are suppressed by the individual "slow-roll" 
parameter, 



1 



2 (HM p f ' 



(32) 



Since (j) 2 is the velocity of a single field, while H 2 is related to the overall density, this quantity is reduced relative 
to its value in the single field case by a factor of N. We can think of ej as the square of the coupling strength of 
each three-point vertex in the interaction Hamiltonian equation (|25[) . The total correction to the power spectrum is 
obtained by inserting in equations (|23|) and (|3"Tj) into (|2~Tj) . 



1 



1 H 



4iV 2 (27r) 3 ej M 2 



2 r tj2 
1 + (ci + 2c2-/Vei)— -j- Ink 



(33) 



Now for slow roll, 



H 

H 2 



N 



2 H 2 M 2 



Nei < 1, 



(34) 



and hence the power spectrum can be rewritten as 

1 H 2 1 

P k 



4(2tt) 3 M 2 e 



H 2 

l + ( Cl + 2c 2 e)— lnfc 



(35) 



which is equivalent to the power spectrum of a single scalar field inflation ip with its one- loop self-correction [42]]. 

The lack of any ./V-dependence in this bound is somewhat surprising. In fact, we will now show that this non- 
appearance is true to all orders in loops, provided slow roll can be assumed. The pair-wise SO(N) symmetry X ! Xj 
in the kinetic term of the action is preserved when we expand the action to higher orders. 8 Consequently, the index 
structure of terms in H mt is always fully pair-wise e.g. <j> I QiQ J Q.j, Q i QiQ j Q.j or <P i 4> j QiQ,;Q k Qk etc. since the 
interactions must preserve this same symmetry. When computing the one-loop diagram, we have to contract through 
the loop to obtain the physical log contribution, in the sense that we have to contract the external lines with at least 
one of the interacting fields in the loop. We have already argued above and in the Appendix [A] that there can be no 
cross-coupling between / and J fields in the four-point interaction, at least at lowest order in slow-roll. Interaction 
terms like 4 l J ( t ) Kd~ 2 (Q I Q I )Q J Q K do exist and can contribute a log divergence, however, these are higher order in 
slow-roll thanks to the extra <j) terms and so their contribution to the one- loop correction will be e/N suppressed. 

Let us now turn to the diagrams generated by the three-point interaction and the even vertex loops. For our 
coherent field argument to be true, the leading 2n-vertex correction must be of the order (Nei) n = e™. We can now 



7 We thank Dan Kabat for a very useful discussion on this point. 

8 The flat target space i.e. G Ir V M 0/ V^j with G IJ diagonal is crucial here. The pair-wise symmetry will not hold if G IJ is not diagonal, 
and we cannot write down the theory as a single equivalent coherent scalar field. Another way of seeing this that one can always redefine 
the fields so that the target space is flat at the cost of generating couplings, both direct and gravitational, in the potential. 
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count the number of diagrams to find the factors of N, remembering that each time a (f) 1 cx Ve^ coupling appears, we 
get a factor of y/l/N from the coherent field relation equation ([34| . Consider the simplest diagram which we have 
calculated in the Appendix: a two-vertex loop with identical interaction term 4> i QiQ j Qj at both vertices. We now 
want to count the factors of N and ej for the correction to the Q 1 propagator. Since we have to contract through 
the loop, we can only contract the Q J of one interaction to the Qj of the other interaction (else we will form a 
disconnected diagram if we contract the J fields at the same space point). Hence, by counting sums, we get a factor 
of N from the J contraction and a factor of ej from the couplings, for a total correction of Nei = e, equivalent to the 
correction for a single coherent field as we have shown in the above. 

Next, consider a more complicated interaction term (j) 1 4> J 4> K Q iQ jQ k ■ Since the fields and couplings have to appear 
pair-wise, there is no 0(<fi 2 ) interaction. There are two ways to contract through the loop, i.e. via Q J and Q , so we 
get two factors of N. However, the three factors of <fi 2 give us the coupling term eiejcx , and equation (|34|) yields 
the correction term A 2 ef — e 3 /N, which is 1/N suppressed compared to the previous case, and hence contributes 
even less. Adding loops will add more factors of e/ into the interaction and thus provide further 1/N suppressions - 
it is easy to see that the extra factors of N coming from the extra loops will never scale more than the suppression 
that comes from e/ = e/N . The point here is clear: since the field labels have to appear pair-wise, we cannot have 
three-point interactions like 4>jQ 2 Q' 1 Q K that could have given us an ~ Ne correction that has no coherent field 
analog. This argument does not depend on the number of loops, since it depends solely on the interaction terms. 

In fact, we can write the N fields in terms of a single scalar field model, where the inflaton is the radial field. 
Assume that each individual field's potential is (l/2)mj(j)j and that they have identical masses mj = m. We can then 
rewrite the fields in polar coordinates, ip 2 = 0i; m which case the Lagrangian becomes 

C= l -{d^Y- l -m 2 ^ + \^{dnf 1 (36) 

where (ip 2 ) ~ (N e/) -1 (H / M p ) 2 . The angular terms thus have large values and damp out quickly, dropping out of the 
inflationary dynamics. Consequently, the set of N coherent fields can be recast as theory with a single scalar field, 
and we are thus unable to derive a bound on N . 

One might worry that the direct coupling term ip 2 (duji) 2 with (dil) 2 — ^(<9f2/) 2 , will generate N—l loop corrections 
that are not scale-free. However, we will sketch below that this interaction can, at the most, generate scale-free loops 
and hence is harmless. Consider the perturbed fields 

V> -> ^ + Q; (37) 

Oj -> Oj + wj. (38) 

Using the symmetry of the interaction, one can show that the leading three-vertex interactions are CljQQuj 1 and 
i/jQuiU! 1 while the leading four- vertex interaction is QQluilu 1 , each with various permutations of derivatives. The first 
three-vertex term can at the most generate a single term, since the attractor solution drives the background angle 
fields to a fixed trajectory and hence we can always pick ilj — (1, 0, 0, 0). For the latter three- vertex term, we have 
to contract through 4 instances of u>i, but it can be shown that the propagator for u>i oc a~ 3 , and hence the loops 
are quickly redshifted away. Meanwhile it is clear from our discussion above that the four-point interaction can at 
most generate scale-free terms. For example, the lowest order diagram is the one-loop one-vertex diagram; since the 
external legs are Q's and the Ui can only contract with itself this is scale free. 9 

IV. DISCUSSION 

In this paper we explore modifications to the dynamics of inflationary models in the presence of N light scalar 
degrees of freedom. We make a distinction between spectator fields, which do not contribute to the background energy 
density, and participator fields, whose potential terms contribute to the inflationary background. As summarized in 
the introduction, a number of very general arguments place finite upper bounds on N, which typically take the form 
N <C M 2 /H 2 . We first consider the dynamics of stochastic inflation in the presence of large number of light fields, 
and show that when N is large there is a distinction between stochastic and eternal inflation which does not apply in 
the single field case. In particular, with a large number of participator fields we show that there is a regime where the 
individual field motion is dominated by quantum fluctuations and thus stochastic, while the overall evolution of the 



9 The analog in the <pi field picture will be the "balloon" diagram. 
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universe is deterministic. Moreover, if the stochastic fields have symmetry breaking potentials, then one can create a 
large number of apparent "pocket" universes, while retaining the ability to end inflation globally, and controlling the 
divergences characteristic of scenarios in which inflation is genuinely eternal. 

Secondly, we explore loop corrections to the 2-point correlation function that provides the inflationary perturbation 
spectrum. Since any light field can run round a loop, these will typically scale with iV. We analyze two subcases - a 
single inflaton with N spectator fields, and N participator fields. In the former case, we find an explicit bound, but 
one which is weaker (by one over a factor that must be small during slow roll) than the simple form described above, 
namely 

M 2 1 

N<^\- (39) 

On the other hand, with N participator fields (A-flation type scenarios), the loop correction is small and independent 
of N. We can understand this result by recasting the action in terms of a composite single scalar field. Finally, in the 
course of this work, we have had need to look closely at the computation of loop corrections in the in-in formalism. 
These calculations raise a number of subtle issues, and we give details of our approach in the Appendices. 

One might whether bounds on TV actually rule out otherwise realistic fundamental theories. Within string theory, 
Vafa [44| argued that T and S dualities ensure that the volume of scalar moduli space is finite, and that there is a 
strict upper bound on the number of matter fields. We are not aware of explicit stringy constructions which would 
saturate the large N bounds described above in any reasonable compactification of string theory and with a finite 
non-zero Newton constant. 

As we have noted, many independent arguments put limits on the allowed value of N. For example, reference [45[ 
notes that if N is too large, the inflaton can decay into a large number of species during reheating, which may cause 
phenomenological problems in the later universe, while [4(| suggests that validity of the perturbative expansion itself 
provides a constraint on N. We have chosen to work in the simplest models of multi-field inflation, where the fields 
have uncoupled potentials. One can multifield hybrid inflation [47] with large number of coupled fields, and it would 
be interesting to ask if such models which are not ruled out by radiative corrections to their potentials but by the 
loop corrections to their power spectrum. Also, while higher correlation functions themselves do not seem to impose 
any bound on N [is} , one could check whether this was also true of their quantum corrections. 
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APPENDIX A: THE "IN-IN" FORMALISM AND ONE LOOP CORRECTIONS TO TWO-POINT 
CORRELATION FUNCTIONS FOR MULTIFIELD MODELS 

In this Appendix, we review the canonical quantization approach to the Schwinger-Keldysh "in-in" formalism 
[H, [2{| and use it to compute the one-loop correction 10 to the 7-th field two-point correlation function. There are 
two types of loops: a two vertex loop of the type considered by Weinberg [29], and a one vertex loop of the type 
considered by Seery [I^] which we calculate below. 

The Schwinger-Keldysh "in-in" formalism [24], [H[ , was first applied to cosmology by Jordan [26| and Calzetta and 
Hu [I?]] • This formalism was reintroduced into the computation of cosmological correlations by Maldacena [28| and 
extended beyond tree-level by Weinberg [29(. We follow Weinberg's notation and methods, with some comments on 
its relationship with the functional method of [26|, [27j ■ The correlation that we want to compute is 

(W(t)) = /(re- , '-^-- H ' tot( * )<B ) t W(t) (Te-^-^ Hi " t{t)dt ^Y (Al) 



See also 49] for a discussion of beyond one-loop effects. 
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where W(i) is some product of fields, Hi nt is the interaction Hamiltonian, both which are constructed out of Heisenberg 
(free) fields and T is the time-ordering symbol. The expectation is taken over the true "in" vacuum. 

One can use the Dyson expansion for the time evolution operator, together with the Baker-Campbell-Hausdorf 
formula to express equation (|Alj) in what may appear to be a more convenient form 29] ; 



(W(t)) = J2* N dtN dt N _ x ... / dt x ([fltatfa), [H int (t 2 ), ... [H int (t N ), W(t)) ...]]) . (A2) 

jy_Q J — oo J— oo J— oo 

While technically equivalent, for these computations this form proves problematic. To perform calculations one 
assumes that at very early times the vacuum state is the bare vacuum, that is, that the interactions disappear. 
Operationally, one implements this by deforming the time contour in equation (|Alj) off the real axis into the lower 
half plane to include a small amount of evolution in imaginary time, killing off the interactions in the far past 
(—00 — > — oo(l + ie)). Unfortunately, in equation (|A2[) this scheme can not be easily implemented. Contours entering 
from the right and the left side of equation (|Aip are treated identically when deriving equation (|A2|1 . However, once 
the vacuum prescription has been specified these contours are in fact complex conjugates of each other, and can no 
longer be freely interchanged. 11 In the rest of this work, we work directly with equation (|A1|) . expanding it to the 
desired order. 

At first order we have, using hermiticity, 

(Q I (t)Q I (t)) 1 =-2Z ( dhiH^it^Q'm'it)), (A3) 



where we have introduced the shorthand oo± = oo(l±ie). Notice that equation (|A3|) is manifestly real. This reality is 
not surprising: we are computing correlation functions and not transition amplitudes as noted in [2H [27j. At second 
order we have, again using hermiticity, 



<Q 7 WQ J W) 2 - -m I dt 2 / di 1 {^ t ) (ti)^ t ) (t 2 )Q J (t)g J W) 

oc_ J — oo_ 

dti / dt 2 (^ t ) (t 1 )Q I (t)g J (i)^ t ) (t 2 )), (A4) 

3 + 



which is also real. Note the time integral contour of the second term. 

At tree-level, the two-point correlation function, {QlQl}, is simply the power spectrum Ql-Ql* ■ This suggests 
that instead of Wick contracting equation. (|A4|) into Feynman propagators, we should contract them into Wightman 
functions instead. Let us see how this works by first defining the contraction 



Q{Qi^QlQi-:QlQi-., (A5) 
where : Q^Qp '■ is the usual normal ordered product and 

J d 3 k Qi - Q J (x,i) - J d 3 ke^ [ aj (k)[/ fe '(i)+aK-kM*(*)] = / d 3 k (Ql+ (t) + Ql k - (t)) . (A6) 



(QkQp) = <[Qk + , Qp"]) = ulu^s'jS 3 ^ + P ). (A?) 



The propagator is now [5l| 



With the contraction, equation (|A5j) . it is straightforward to prove that the correlations of equations (|A4j) and (|A3j) 
are a sum of all possible contractions into both connected and disconnected pieces as per the usual Wick's Theorem 
in standard quantum field theory. We ignore the disconnected pieces, i.e. the "vacuum fluctuation" pieces, in which 
vertices are connected only to other vertices and no external lines. Since we are computing correlation functions 12 , 
the pieces automatically cancel [29l |. 



11 Note that this issue is not manifest at tree level, and the conclusions of [29l |34H are robust, other than with respect to the coefficient 
of the N spectator loop correction, as discussed in Section III 1 1 Musso |50| has developed a diagrammatic formalism for correlation 
functions in the "in-in" formalism. However, since this is based primarily on equation jA2t . one would need to check it carefully before 
employing it for an explicit calculation. 

12 Whereas when we compute transition amplitudes, they contribute an overall phase 1531 . 
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As an aside, we note here that in original "in- in" formalism of Schwinger-Keldysh (and see also [26|, [U [HI ) , the 
original fields are split into + forward time fields and — backward time fields each with their own generating functional. 
One can then compute all four possible Green's functions for the fields (+, +), (— , — ), (+, — ), (— , +), and then all the 
possible contractions of the correlation will be one of the four above. In other words, the doubling of the fields into 
+ and — sets provides a convenient book-keeping method for keeping track of the contours. In our formalism, we 
have a single contraction equation (|A5j) , but we pay for this simplicity by with having to explicitly keep track of the 
contours of the integrals we must perform. 

Our strategy is as follows: 

• Expand equation (|A1[) to the desired order keeping careful track of the vacuum prescription contours. 

• Insert the interaction Hamiltonian H lnt and expand all fields using equation (|28D . 

• Expand using Wick's theorem and the contraction defined in equation (|A5|) discarding disconnected diagrams. 

• Perform the integral over time(s) leaving only integrals over the internal momenta. 

• Regularize the remaining integrals to obtain the final answer. 

One should note that these loops contain both UV and IR divergences, unlike the case of [2j|. The presence of the 
IR divergences means that our use of dimensional regularization will yield incorrect finite terms. However, we are 
only interested in the logq dependence, which is correctly computed by dimensional regularization. In this paper, we 
are concerned with large- TV effects, but a clear understanding of the mechanism which regulates these logs is clearly 
of critical importance, and we plan to pursue this topic in a future publication. 13 



1. Two- vertex loop 



The two-vertex loop is generated by a three-point interaction term [3S 



—V^QiQjQj + - 1 rV2^id- 2 QiQjd 2 Q u 



(A8) 



with the coupling term, 



2H 2 ' 



(A9) 



where here, and in all subsequent calculation in this appendix, we have dropped all factors of M p to simplify notation. 
We have also dropped the sums over the J and I since all the fields are identical - most of the diagrams have identical 
amplitudes and we will sum them later. 

For the two- vertex terms generated by equation (|A8|) . there are two types of vertices, and we have to compute the 
contributions from all the possible combinations of two vertices and internal loops (see Fig. ([1]), giving us a total of 
six diagrams. In the following, we sketch the derivation for the first diagram of Fig. ((T|) with QiQjQj vertices at 
both ends - the other diagrams are computed similarly. Using the the first term of equation (|A8|) in equation (|A4[) . 
and Wick expanding everything with the contraction of equation (|A5[) . we get 



a 2 (n ) v/277(n) a 2 (r 2 ) T Tl 



3(2?r) 9 5> 



dr 2 



dri- 



J d 3 k J d 3 fc' / f7 fe J (T 1 )f7/*(r 2 )^(r 1 ){7 fc J ,*(r 2 )5 3 (q+k+k') 



One can in principle impose a horizon cut-off, as suggested by Lyth in Ref. [53 | and applied in [55l |. Boyanovsky et al. [56t 157. 58, 59] 
have suggested that the IR divergences are regulated by the slow roll limit. However, their approach requires one to analytically continue 
a combination of the slow roll parameters, which are physical, and in principle measurable, quantities. 
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dri 



dri- 



er 



! (r 1 )^/2i7Ra 2 (T 2 ) v /2i7R rr/ 



4 4 

'3iJ f TTJ f~ \ttJ*{~ \t'tJ*(~- \t'tJ*(~ \ ^ 3 / 



^(Ti)^*(r)f/^(r 2 )^(r) 



x / d 3 fc / d 3 fc"C/ fc J (r 1 )C/ fe J *(r 2 )C/^(r 1 )^*(r 2 ) ( 5 3 (q + k + k') 



(A10) 



where we have dropped all disconnected terms. An overdot in this expression denotes a derivative with respect to 
conformal time, r. 

We now integrate over the times first, since the interactions of equation (|A8|) satisfy the late time convergence 
conditions 14 . In near de Sitter space H « const, and the mode functions Uf can be approximated by 



Ul{r) 



H 2 



2(27r) 3 fc 3 



(1 + ik T )er lkT . 



(All) 



Although in general the couplings ej(r) are not constant, since they are changing only slowly we make the further 
simplifying assumption that they are roughly constant at late times. Plugging equation (|A11[) into equation (|A10p . 
and then integrating over the times, we obtain 



J d 3 x e lq(x - x ' ) (vac,in|Q / (x,T)g / (x,r)|vac,in 
-^3^2 J d 3 k d 3 k'S 3 (q + k + k') jej 



16(2 
+ 0(qT)+ ...} 



5 kk 3 kk 



4 q 7 K 4 q 6 K 2 



kk'{K + q) 2 



(A12) 



where we have assumed identical fields hence ej = ej and used K = k + k' + q. 

Considering only scales that are larger than the Hubble horizon gr < 1 allows us to also drop the second to last 
term in equation (|A12[) . The rest of the 2-vertex diagrams can be computed identically, and the final answer we 
obtain is 



J d 3 xe lq (x - x ' ) (vac,in|Q / (x,r)g / (x',T)|vac,in) 2 
J^ N Jd 3 kJd 3 k'S(k + k' + c l ) 



39 kk' 5 k' 



1 

K 2 

1 

IX 1 



51 kk' 31 k' 5 k' 3 k' 



4 q 6 2 k 3 q 2 8q 4 k Aq^k 2 



3 k 1 kk k 
2 kq 2 ~ 2^ ~ 2 Wq 



1 



kk' 



31 k' 


+ 4 q 3 k 3 _ 


2k'' 




qk 3 





(A13) 



Our task now is to compute the physical logs from equation (|A13|) . From simple dimensional analysis, they take 
the form 



mull 



d 3 kd 3 k' (5 3 (q + k + k') 



k m k 



3+/+m— n+6 



(q + k + k') n 

To obtain the coefficient Fn t7n>n \, we use the identity 



F(i,m,n) m 9 + polynomial divergences. (A14) 



/rco pk+q 
d 3 kd 3 k' 6 3 {q + k + k')f(q, k, k') = 2tt / kdk \ k'dk' f(q, k, k') 
Jo J\k-q\ 



(A15) 



We can extract the coefficients of physical logs from the above expression using l'Hopital's rule. We find the limit for 
both terms going to infinity by differentiating the right hand side of the above identity and equation (IA14[) as many 



14 See also [6011 . If we switch the order of integration, we end up swapping an ultraviolet divergence in the integrals over internal momenta 
for a divergence in r. 
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times as needed to tease out the log divergences, and comparing the result to the same operation applied to the right 
hand side of equation (|A14[> . 

Dropping all other polynomially divergent terms and the log IR divergences, each term in equation (|Al3j) contribute 
the following: 



Fi 



(1,1,1) 



—TT TP 

~ "16"' t ( 



1,1,2) 



^(1,-3,2) — 7r , ^(1,-3,3) 



2tt, 



■£(1,-1,4) =0, £(1,-2,2) = 7T> ^(1,1,4) = f ) £(1,1,3) = T"> 
£(1,-2,3) =0, £(1,-1,1) = -jr, £(1,-2,1) =71", -£(1,-3,1) =0 ' 
£(1,-1,2) = TT- 

Given identical fields, such that ej = ej, then putting together everything give us the total contribution from the 
two-vertex one-loop diagram: 



J (fx e iq ( x ^ x ' ) (vac,in|Q / (x,T)Q / (x, r)|vac,in) 2 
where denotes scale-free polynomial divergences. 



H A Ntj 
2(2vr)V 



( ) 7r In q 

V 120 J 



(A16) 



2. One- vertex loop 

On the other hand, the one- vertex loop is generated by a four-point interaction term derived in Appendix [B] 

1 



/ d 3 xa 3 Y< 

i,J 



iH ^d i Qjd i Q J d- 2 (d J Q I d :l Q I + Qid 2 Qi) 



where, 



+-^QjQjd- 2 (d l Q I d J Q I + Qid 2 Qi) 

4/7 

+-^d- 2 (d j Qjd J Qj + Qjd 2 Q J )d- 2 {d, J Q I d ] Q I + Q z 9 2 Qj) 

4/7 



^f3 2 , 3 ^ d- 4 (djdkQ'dkQj + d 3 Q I d 2 Q I - 9 2 Q J ^Qj - d^d^Qi 



(A17) 



(A18) 



The four-point interaction (|A17|) is explicitly independent of the background potential 15 V . This means that the 
one-vertex one-loop correction to a field / from all other fields J ^ I is the same for whether or not J are spectators 
or participating fields. Each four-point term has the form Q 1 Q 1 Q J Q J , meaning that depending on the external lines, 
we can contract it with either the J or / fields, thus each term will effectively generate two different diagrams. 

Fortunately, only the self-interaction term, i.e. 1 = 3 contracted with / external lines, contributes a physical log. In 
the UV, all other terms diverge polynomially and we assume that they can be absorbed by renormalization. Heuristi- 
cally, this is because the interactions are secretly mediated by gravitons and this, combined with the symmetry of the 
action, prevents any non-self-interaction loops from contributing, as described in Section IIIII From a diagrammatic 
perspective, to yield a log divergence, the final momenta integrand must possess a scale. That is, it must have the 
form ~ fc t 7(k±q)' 3 for some a £ Q and (3 £ Z + , where q and k are the external and internal momenta respectively. 
However, note that the fields operated by d~ 2 is always identically paired i.e. they appear as d~ 2 (Q I Q I ) and never 
d~ 2 (Q I Q J ). Hence if I ^ J, the integrand can only possess momentum factors like l/(k 2 ) or l/(q + q') 2 , which only 
contribute polynomial divergences. For example, a four-point term with an interaction (dropping time derivatives as 
they do not affect the final result) Qjd~ 2 (QiQi)Qj for a (Q I Q I ) correlator will yield the diagrams shown in FigJ5]if 
I =/= J, which is simply a vacuum fluctuation diagram multiplied by a propagator. 

We now turn our attention to the self-interaction term, where I = J interaction terms are contracted with / external 
lines. This calculation is operationally the same as that done by Seery [42j for the single scalar field case using different 



The dynamics of the field depends on the potentials and their mode functions, i.e. their green's functions will differ, but the coupling 
terms have the same structure. 
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{Qiix)Q I (x , )d- 2 (Q I (z)Q I (z))Qj(z)Qj(z)) oc 



(A19) 



+ 



FIG. 2: The non-self-interaction four-point loops factor into a vacuum fluctuation piece times a propagator. The ellipses 
indicate contractions which lead to polynomial divergences. 

techniques. Fourier transforming equation (|A17|) and switching to conformal time, the interaction Hamiltonian 
becomes 



(4) 



J -°° I.J 



a °"( k :P) r\l r\l r 



o-( k .P) AI^I, 



AH (k + p) 



kQp<#Qb*(k + P + a + b) - — a • b ^^ 2 QiQ T p QiQi6(k + p + a + b) 



2 ( 3 o-(k,p) o-(a,b) z(k,p) ■ z(a,b) o z(k,p) ■ b^ Al „ lAJ 



+ 



QkOpQa 7 Ob^(k + P + a + b) 



v 4(p + k) 2 (a + b) 2 (k + p) 4 (a + b) 2 (k + p) 4 

where overdot again denotes derivative with respect to conformal time. We use the notation of [42| . 

z(k,p) = cr(k,p)k-cr(p,k)p, 
cr(k,p) = kp + pp. 



(A20) 



(A21) 
(A22) 



After some tedious but straightforward calculation reminiscent of the previous section, and considering modes well 
outside the horizon, qr — > 0, we find that the one-vertex self-correction is 



(Q I (t)Q i (t)) 1 = (2vr) 3 / d 3 k[A s +B s +C s ], 



with the following contributions 



A, 



H 4 



32(271)^ 



y q - k 



Qq 2 5\ fc 2 + q-k 10 q 2 + k • q 



k 3 k) (k + q) 2 k (k + q) 2 



(A23) 



(A24) 



and 



where 



B„ = - 



H 4 



32(2tt)V 



1 k 2 + q • k 5k q 2 + k • q 
k (k + q) 2 + 7 (k + q) 2 



C.s = 



tf 4 



5k (2a 2 — k 2 ) fi 

/(k, q, -k, -q) + 1 p j /(q, k, -q, -k) + -/(k, q, -q, -k) 



/(k,p,a, b) 



8(2tt)V 

3 er(k, p) cr(a, b) z(k, p) • z(a, b) z(k, p) • b z(a, b) • p 



(A25) 



(A26) 



<5(k + p + a + b). (A27) 



^4(p + k) 2 (a + b) 2 (k + p) 4 (a + b) 2 (k + p) 4 (k + p) 4 
The integral in equation (|A23j) is divergent and so it needs to be regularized. After regularization we obtain 



H 4 5tt 
2(27r) 3 q 3 16 



In q , B s 



H 4 TT 
~2(27r) 3 g 3 48 



lnq , C s 



H 4 TT 



2(27r) 3 g 3 3 



ln 9) 



(A28) 
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yielding the following log q contribution 

{Q I {r)Q I (r))=-^^\nq + ..., (A29) 

which is our final answer. 

Finally, we sketch the technique we used to extract the logg terms from equation (|A23[) . The key idea is to write 
the integrals like 



d k p«'(k±q)2/3' (A30) 
(where a is a half integer and (3 £ {1, 2, 3}) as a sum of terms of the form 

/^)/^ V-n(kU^ - (A31) 

We can then use standard techniques (e.g. f52}) to evaluate the integrals, in combination with the following trick. 
Define 

^ (^^V (A32) 
(l)a,/J = (A33) 
Now since k l is integrated out, the only vector quantity left is q l so the following must be true 

= B^pq* , (fc*fc J ') ,/3 = C^^'+^W (A34) 

where B a< p, C a .p and Z? Qj| a are coefficients that may contain ultra-violet divergent components. We can then dot 
both sides of equation (|A34[) with g's and complete the square to eliminate the numerator (q • k) terms. For example 
(for the B term), 

B a3 q 2 = ±i<(q ± k) 2 ) Qi/3 T ^ + <7 2 ) Qj/3 , (A35) 

and the first term cancels one of the powers of (k ± q) while the second term cancels out one power of k 2 . This leaves 
us with 



B a ,0 — ±2^-fa,/3-i T 2^a,f} T 2~~2-k>-i.0j (A36) 



We can iterate this trick until we remove all terms with powers of (q ■ k) in the numerator or cancel all the (k ± q) 
terms in the denominator, resulting in integrals that are polynomially divergent and hence can be discarded. The 
remaining integrals are in the form of equation (|A31|) and can be easily regularized. The C a $ and D a ^ terms can 
be similarly computed by dotting twice with q's. 



APPENDIX B: MULTIFIELD 3-POINT AND 4-POINT ACTION 

In this paper, we make use of the Arnowitt-Deser-Misner (ADM) formalism [6l[ to expand the action equation (JTJ 
to 4th order in perturbations. The derivation is straightforward if rather tedious (see for example refs. [33|, E3] for a 
detailed application of this formalism), so we simply collect the results. The background N- field action is 



5' 



d 4 Xy/g 



R 
~2 



E 



i 



(c%) 2 + vH0/ 



(Bl) 



the ADM metric is 



ds 2 



-N 2 dt 2 



iij(dx l 



N i dt){dx j + N j dt), 



(B2) 



and we choose to work in the spatially flat gauge so hij = a 2 (t)6ij. In other words our metric perturbation has been 
set to zero by a gauge choice. In addition, we focus on the scalar perturbations, and ignore vector and tensor pieces. 
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This means that the diagrams we computed do not have graviton propagators or loops. In this gauge the fields have 
non-zero perturbation 



Qi 



(B3) 



Field indices are summed over when contracted, X I Yi = J2i XiYi and V = J2i Vi{4>i)i with no cross-coupling terms 
between the fields. 

The quadratic action is (exactly) [33| 



S> = ~ / did 3 " <■ ■"• 



Qi - "-'u)(J l r-[Vjj-~(^ I h 



lr>J 



Q Q 



The third order action is, to leading order in slow roll [331 ], 

1 



dtd 3 xa 3 
1 SL 



a 3 SQ 



where the 6L/SQ J is the first order equation of motion. 

Finally, the fourth order action is, to leading order in slow roll, 



Si = I dtd 3 xa 3 



j (V 2 (d j Q I d j Q I + Q I d 2 Q I 



- ~&,id 2 /?2j + X2(d l Q I d l Q I + Q I 2 Q I ) ~ Q 1 ' (i^Qi 
where the auxiliary fields \i and 02,j a- r G, to leading order, 

iftj ~ d- 4 (djdk&dkQi + d 3 Q I d 2 Q I - dP&djQi - d^djdmQ! 



(B4) 



(B5) 



(B6) 



(B7) 



d 2 X2 = -^ d >Q' d >Qi - \ d ~ 2 (djQ'djQi + Q I d 2 Q 1 ) - -^Q'Qi- 



(B8) 



Note that the 4-point action equation (|F36|) will generate N 2 non-trivial diagrams, since there are two sums over the 
field indices. Note also here that this expression reduces to those of [42j in the single field limit. 

APPENDIX C: QUANTIZATION OF THEORIES WITH DERIVATIVE INTERACTIONS 

In this appendix, we describe the procedure we use to canonically quantize the classical theory, as we encounter 
Lagrangians with interactions containing time-derivatives of the fields. A pat h integral formalism is given in [62j ■ For 
theories with up to 2nd order in time derivatives, a treatment is given in 16311 . In the problem we are considering, we 
encounter interactions up to 3rd order in time derivatives, so we extend [63[ at least to the next to leading order in 
slow-roll. Extension to all orders is straightforward which we will leave for future work. A path-integral approach is 
pursued by Seery [42| . 

We follow the usual procedure in canonically quantizing a classical theory specified by a Lagrangian density, C(Q, Q). 
That is, we define the momenta conjugate to the field Q by; 



- 5£ 

and construct the Hamiltonian density, 7i, as the Legendre transform of the Lagrangian density; 



(CI) 



n = nQ-c, 



(C2) 
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where Q is expressed in terms of it. We then move to an interaction picture by separating the Hamiltonian into its 
quadratic part Hq and higher order part 7ii nt and replace ir in 7i; nt with the interaction picture 717 given by 



Q = 717 



d-K 



(C3) 



The question we want to address here is, what is 7ii n t? Naively, one might guess that 7i; n t = — Ant ; where 
Ant = £ — £q and £0 is the quadratic part of C If the only time derivatives of the field are in a canonical kinetic 
term, this is certainly the case. However, when time derivatives are present in the interaction terms, these can modify 
the relation between ir and Q, and the construction of the Hamiltonian then generates extra interactions. Fortunately, 
at the order we are working the additional terms generated are either subleading in slow roll, or higher order in the 
fluctuations. To see this, note that the Lagrangians we consider have the schematic form 

C =\Q 2 - V(Q) + (sfeh + / 3 ) Q + \ {Vegi + .92) Q 2 + ^h,Q 3 + 0(eQ 3 ) + 0(eQ 4 ) + 0(Q 5 ) (C4) 

where e is the usual slow roll epsilon and the subscripts of f m and g m denote a term of order m in fluctuations, Q. 
The terms containing no time derivatives are gathered into V(Q). To proceed, we assume that \Q\ ~ \Q\ ~ |tt|. A 
straightforward calculation then shows that 

n = ^ + V(Q) - {V~ef2 + h) Q-\ {V~egi + 32) Q 2 - ^Q 3 + e (j 2 giQ + \glQl 2 ) (C5) 

+0(eQ 3 ) + 0(eQ 4 ) + 0{Q 5 ). 
= H -d + 0(e Q 3 ) + 0(e Q 4 ) + 0{Q b ). (C6) 

So, to leading order in slow roll and to fourth order in fluctuations, it is safe to take TL m t = — Ant) the correction 
being at the most of 0(e) and thus subleading during inflation. 
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